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e-j . Abstract 

On" We suggest a parameterization of the matrix element for B° — > D* + ujit decay 

Q_i' using kinematic variables convenient for experimental analysis. The contributions 

of intermediate urn- and D**-states up to spin 3 have been taken into account. The 

angular distributions for each discussed hypothesis have been obtained and analysed 

_ -, ' using Monte Carlo simulation. 
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1 Introduction 

oo ■ 

The discovery of excited -D-states (referred to as _D**-states) stimulates interest in their 
spectroscopy and D** — > D^n decay properties. There are four P-wave states, which 
are usually labeled D* (j£ = 0+ /2 ), D[ (Jf q = 1+,), D 1 (j£ = l+ /2 ), D\ (Jf q = 2+ /2 ), 
where J is the spin of the meson and j q is the total angular momentum of a light quark 
q = (u, d), which is the sum of the orbital momentum / and the light quark spin s q . In the 
heavy quark limit, the angular momentum j q is a good quantum number. Conservation 
of parity and angular momentum imposes constraints on the strong decays of the D** to 
D^n. Two states with j q = 1/2 decay to the D^Tr-st&te in S-wave while two other with 
j q = 3/2 decay in D-wave. Since the decay width T ~ Q 2i+1 , where Q is the magnitude 
of the daughter particle momentum, L is the orbital momentum between decay products, 
and Q is small, Di and D\ have small decay width of about 20 MeV, but Dq and D[ are 
expected to be quite broad with decay width of about 300 MeV [TJ |2]. 

A further study of these states will allow a more in-depth comparison to be made 
with theoretical predictions such as Heavy Quark Effective Theory (HQET) (3j H] and 
QCD sum rules |5]. The last experimental studies of D** mesons were performed in 
B~ £>M + 7r~7r~ [7] and B° D ir + 7i- [8] decays. These states have also been 
studied in semileptonic S-decays [9] . Thus, understanding of their properties is significant 
for reducing uncertainties in the measurement of semileptonic decays and determination 
of the Cabibbo-Kobayashi-Maskawa (CKM) [TD] matrix elements \V c b\ and \V u b\. 

-D**-states can be also produced in other hadronic B decays, e.g., B — > D*utt. Here, 
D** production is described by the W vertex instead of the transition Isgur-Wise func- 
tions j3], which describe these states in the D^nir modes. This channel was first observed 
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by the CLEO [TT] and BaBar [12] collaborations, the latter finding an enhancement in 
D*7i mass due to the broad _Di(2430)°-state, representing a P-wave of a D meson. 

Let us note that light mesons decaying to the con final state (e.g., p(1450), 61 (1235) 
and their excitations) appear in the color-favored mode of this process. Thus, a possible 
contribution of these resonant structures to the total branching fraction can be measured. 
The p(1450)-resonance, dominant in this mode, was observed by both collaborations [HI 
IT2] . but the 6i(1235)-state was not observed in this channel. 

An amplitude of three-body decay can be written as a sum of the contributions cor- 
responding to the quasi- two-body resonances [El El |8] . Analysing experimental data one 
has to determine relative amplitudes and phases of different intermediate states. To do 
this, one needs the amplitudes expressed via kinematic variables convenient for Dalitz 
plot analysis^. These expressions can be used for optimization of selection criteria and 
creation of efficient Monte-Carlo generators. 



2 The general method 

A weak P(0~) — y R(J P )1~ decay amplitude (for J > 0) includes three independent terms 
while a strong R{J P ) — > 1~0~ amplitude can have one or two independent terms. We 
can parameterize a decay matrix element using a set of different independent bases. In 
general, we can use the basis of covariant amplitudes or helicity basis etc. Since the real 
particles Di and D[ are expected to be close to the pure j q = 3/2- and j q = 1/2-states 
and their decays have particular orbital momenta, it is convenient to use the basis of 
amplitudes describing decay with fixed angular orbital momenta in the B and resonance 
rest frames. 

In this paper we use an isobar model formulation in which our decay is described 
by a coherent sum of a number of quasi-two-body amplitudes. The amplitudes can be 
subdivided into two channels. The effective Hamiltonian for Cabibbo-favored decays can 
be reduced to the color-favored and color-suppressed forms [131 E] : 

Hcf = ^V^McT^iuT^ + ^H'J, 

Hcs = ^7^(05,(^^6) + dflS), (1) 

where Gf is the Fermi constant, C\ and C2 are the Wilson coefficients and T M = 7^(1—75). 
The coefficients a\ = C\ + C2/N and a 2 = C2 + C\/N, where N is an effective number 
of colors. The terms Hi = \ j^ a=l {c\ a Y ^b)(u\ a Y^d) and Hi = \ ELi^AT^XmAT^) 
(A a are the Gell-Mann matrices), involving color-octet currents, generate non-factorized 
contributions. The other non-factorization source is the non-factorized matrix element 
of the product of the color-singlet currents. It includes loop current- current terms. The 
color-favored and color-suppressed channels are shown in Fig. [TJ We show tree diagrams 
only, however, not all the intermediate states are described by them. In this paper we do 
not apply the factorization method but consider all intermediate resonant contributions 
up to spin 3 allowed by the momentum-parity conservation. The color-favored term 

1 In the case of decays with more than three particles in the final state, the term Dalitz plot is used in 
a general sense to refer to the distribution of the chosen degrees of freedom used to describe the decay. 
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Figure 1. a) Color- favored and b) color-suppressed channel. 



receives a contribution from the o;7r-resonances, e.g., p(1450) and &i(1235). Since these 
resonances are broad, this channel allows factorization to be precisely tested [15]. The 
color-suppressed term receives a contribution from the P**-states, which are P- and D- 
wave excitations of the cu states. 

Let us consider briefly the spectroscopy of the P-wave cu excitations. We have Jj u = 
I3/2"' J-H = ^3/2"' Jfv. = ^5/2"' an< ^ = ^5/2 states. Again, as discussed above, two states 
with j u = 3/2 decay to the D^ir-st&te in P-wave and two other with j u = 5/2 decay in 
F-wave. 

Observable cu-states with the same J p = 1 + (J p = 2~) quantum numbers are two 
linear combinations of pure j u = 1/2 (j u = 3/2)- and j u = 3/2 (j u = 5/2)-states. Thus, 
the physical D\ and D^-states are as follows: 

\D X > = sin^i \j u = 1/2 > + cos^i e~ m \j u = 3/2 > , 
\D[ > = cos^x \j u = 1/2 > - sin^x \ Ju = 3/2 > , 

where $1 and $2 are mixing angles. 

Let us discuss kinematic properties of the considered process. In the final state we 
have six particles, namely, D° and tt + from the D* + decay, 7r + , 7r~ and 7r° from the cu 
decay and ix~ from the B° decay. The B° decay is described by two invariant masses 
squared of the D*n ijn 2 DtlT ) and urn (m^) systems, the one corresponding to resonance 
mass labeled as q 2 . 

The u decay is described by five variables. We use invariant masses squared Mq = 
(P+ + P_) 2 and M\ = (P + + Pq) 2 (here Pj is a 4-momentum of the pion ir l from the cu 
decay, % = ±, j§, the azimuthal angle of the 7r° in the u decay plane, and two angles 
(polar 9 and azimuthal 0) for a vector n normal to the u decay plane. Let us note 
that the u> — > tttttt decay proceeds through two mechanisms. The first one involves an 
intermediate p-meson. Experimental studies of the e + e~ — > 3tt reaction have confirmed 
the Gell-Mann- Sharp- Wagner suggestion [16] that the u — > 3n transition is dominated by 
this contribution. The second mechanism represents the non-resonant contribution. This 

2 The u! invariant mass squared is p 2 = (P + + P_ + Pq) 2 - 
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contact contribution can not be excluded because interference between these mechanisms 
leads to a sizeable effect in the decay rate. 

The D* + decay is described by two variables. We use polar and azimuthal ip angles 
for the D° momentum in the D* + rest frame. For further applications we assume the width 
of the D*-meson to be negligible (T D *+ 0.1 MeV <C m D *+ — (m D o + m w +) rs lOMeV). 
To describe the intermediate resonance decay, we use polar £ and azimuthal ( angles 
for the daughter particle momentum in the resonance rest frame. The polar angle £ is 
expressed via the invariant mass squared m 2 Dt7T for the cu7r-states and m 2 ^ for the D*n- 
states. Moreover, the matrix element does not depend on the azimuthal angle ( for the 
unr- as well as for the _D*7r-states. 

A further definition of angles depends on the decay channel. Figure [2] shows the decay 
scheme and definition of the angles for the o;7r-resonances. 




Figure 2. Complete visual definition of the angles for the ow-resonances. The angles 6 and (j) 
are defined in the uj rest frame, the angles j3 and ip are defined in the D* rest frame and the 
angle £ is defined in the urn rest frame. 

Figures |3] and H] define these angles using momentum variables for the uir- and D*7r- 
resonances, respectively. The notations are as follows: the variables p, Q, I, q are the 
four-momenta of the u-, D*-, D-meson and an intermediate resonance, respectively, while 
p, Q, 1, q are the magnitudes of their three-momenta in the mother particle rest frames. 
In Figs. [3l H] the directions of these momenta define angular variables 6 and <fi in the u 
rest frame, /3 and ijj in the D* rest frame and £ in the resonance rest frame. 

In this paper each compound particle is described by a relativistic Breit-Wigner (BW) 
with a g 2 -dependent width. Such an approach is not exact since it does not take into ac- 
count final state interactions and is neither analytic nor unitary. Nevertheless, it describes 
the main features of the amplitude behaviour and allows one to find and distinguish the 
contributions of different quasi-two-body intermediate states. Thus, the denominator of 
the BW propagator is: 

D R (q 2 ) = q 2 -m 2 R + im R T R (q 2 ). (2) 

It corresponds to the intermediate resonance R with mass m R and g 2 -dependent width 
Fr. The numerator of the propagator is to be the sum over polarizations of the resonance 
and depends on its spin. 
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Figure 3. Definition of the angles for the w7r-resonances. Color-favored channel, a) The ujtt rest 
frame, b) the D* rest frame and c) the u rest frame. 




Figure 4. Definition of the angles for the D** -resonances. Color-suppressed channel, a) The 
D** rest frame, b) the lv rest frame and c) the D* rest frame. 



6 u;7r-resonances 

We consider such o;7r-states, which can be combined to J p = 0~-, J p = 1 + (6i (1235))-, 
J p = 1- ( p (1450))-, J p = 2"-, J p = 2+- and J p = 3-(p 3 (1690))-states. Let us note that 
J p = 0~, J p = 2~ and J p = 2 + charged states, which decay to the ui7r- final system, have 
not yet been observed at the present time Such states have the isotopic quantum 
numbers I = 1 + . It is natural to assume that these states are members of the b and 
p-families. 

The matrix element for production of the J p = 0~ intermediate state (labeled as po) 
is given by: 

G p 



M 



B^D*p 



V2 



V cb V: d < D*p \(uT ll d)(cT»b)\B > 



(3) 
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Parameterizing this amplitude in the covariant form, we haveJl 



M B ->d* P0 = ^V cb V: d g po F P (q 2 )(e*q), (4) 

where e p is a polarization vector of D*, g Po = ai/ po 0, f Po is a weak decay constant of the 
Po and Fp(q 2 ) is a transition form factor. 

The strong amplitude for the pcrdecay is presented as follows: 

M po ^ wn = 9 P0 ^Fp(q 2 J p 2 )(v*q) J (5) 

where v p is a polarization vector of u, g pou)n is a coupling constant and F P (q 2 ,p 2 ) is a 
transition form factor. The amplitude describing the u decay comprises the contributions 
from the intermediate p- meson and 3n phase space: 

M^ 3n = 9ujpn (p 2 ) (a 37r + D ^m 2 )Z[M 2 ) J V A (P, p +, p o)( nv ), ( 6 ) 



i=±,0 



where 



(7) 



v/A(p,P + ,P ) 

is a unit 4- vector normal to the u decay plane and A(p, P + ,Pq) is the Kibble determinant. 
Other notations are described in the Appendix. Here and further e pupa is the Levi-Civita 
symbol and e i23 = +1- The amplitude corresponding to the D* decay is 



Md*^Dtt — 9D*Dic{£l)- 



The factor 



g D * D «(Q W(p ) ^ + ^ D p i(M 2 )z(M 2 ) ) d d 4Q 2 )d u (p 2 ) 1 (9) 

is common for all intermediate states, and w-decay part can be expressed via the phase 
integral W(p 2 ), presented in the Appendix. 

The total rate for B — > D*ujtx decay expressed via the branching fraction Bd*+~^d°-k+ 
and the phase integral W{p 2 ) can be presented as follows: 

dT = ia mo 2 r^T~ m / 2M2 d P i d cos 6d(f)) (dcospdip) (dq dcos£), (10) 



Q 

where ms is a I?-meson mass, and the matrix element M describes particular dependencies 
for the different intermediate channels. 



3 Here the term (e*Q) is neglected because the longitudinal currents arise far from the resonance, where 
they should be suppressed by transition form factor behavior. However, they also modify the angular 
dependence of the amplitude. Throughout this paper the longitudinal currents are neglected. 

4 The coefficient a\ is expressed via Wilson coefficients, as discussed in the previous section. 
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The matrix element for the B — > D*Rj transition, where Rj is the intermediate 
resonance with the integer total spin J@, can be parameterized in terms of the amplitudes 
with the definite angular orbital momentum L as follows: 

GVt " ' \cje^e'< J h:q p Q a F L= j(q 2 ) + 



M, 



B^D'R, 



/^V cb V ud gj 



(e'^)Q)^q))F L= j. 1 ( q 2 ) + 



fj,j-i(q 2 j 

+ iC J+1 ((e'*WQ)(e*q) - fj,j + ,{q 2 ){e^e*))F L=J+1 {q 2 ) 



Here, Cj is the relative amplitude, which is in general complex; F^q 2 ) is a transition 
form factor corresponding to the orbital momentum L\ g% = a\f R , f R is a weak decay 
constant of the vector resonance, when q 2 = m R , g 2 = m BdBD* R2 an d 93 = ^b9bd*r 3 are 
appropriate coupling constants; e ^ is a convolution of the resonant polarization tensor 
of rank J and momentum (j^ 



e '(J=i) = £ ' - 



' (J=2) = £ 'uaQ a /' m B, < (J=3) 



/j,J±l(<? 2 



2 r»2 



m B _ m D* - I 2 + 2aj,j±i77lD* VT 



(12) 
(13) 



and 



a 1(Q = -1, a 1)2 = +2, a 2 ,i = -1, ^2,3 = +3/2, a 3j2 = -1, a 3j4 = +4/3. (14) 

The parameterization of the matrix element describing the resonance decay depends 
on its J p quantum numbers. Thus, resonances with J p = 1~, 2 + , 3~ are described by the 
following matrix element: 

M Rj ^ n = ~gje^e'yh;q p p a F L=J (q 2 , P 2 ), (15) 

where g x = g RlU7n g 2 = m R g R2UJ7T and g 3 = m 2 R g RsL07r are appropriate coupling constants; 
Fl{q 2 ,P 2 ) is a transition form factor and 

et l) = ^ ^'t 2) = V a /rn R , e'^ = s>V/m|. (16) 
The discussed resonances with J p = 1 + , 2" are described by the following matrix element: 



M. 



R j— turn 



93 



Cj^mU(e (J) v*) 



fj,j-i{q 2 ) 



.e^p)(v*q))F L= j^(q 2 , P 2 )+ 



+ C J+1 ((e'^p)(v*q) - f w {q 2 W^v*))F L = J+ i{q 2 ,p 2 ) 

where Cj is the relative amplitude, which is in general complex, and 

2q 2 p 2 



fj,j±i(q 



q 2 +p 2 - m 2 + 2aj j j±i-\/p 2 q 2 ' 



(17) 



18) 



where m is the charged pion mass. 

Then we move from the covariant amplitudes to the expressions depending on the 
selected angles, which are defined in the intermediate particle rest frames. 



5 As emphasized above, in this paper we discuss resonances with J = 1, 2, 3. 
6 The notation e ^ is not related to the resonance helicity state. 



7 



4 J D**-resonances 



The decay rate for the channel with P**-resonance production has a form similar to ( JTUl) . 
As already mentioned, in this case the angles (9, 0, £, 0, ip) differ from their analogues for 
the utt states and are described in Fig. HI Here we discuss two J p = l + -states and a 
J p = 2 + -state, which correspond to P-wave in the spectroscopy of the cu excitations as 
well as a J p = l~-state, two J p = 2~-states and a J p = 3~-state corresponding to the 
D-wave excited cw-states. Pure J P u = l+, 2 (J P u = 2g, 2 )- and J P u = lt, 2 (J P u = 2^ 2 )-states 
decay to the D*tt in S- (P-) wave and D- (F-) wave, respectively. As discussed above, 
observable J p = 1 + (J p = 2~) states can be a mixture of pure j u = 1/2 (j u = 3/2) and 
ju = 3/2 (j u = 5/2) states. This fact has to be taken into account for the total amplitude 
construction. The parameterization of the matrix elements for all P**-states is similar to 
the case of the cu7r-states. However, mutual substitutions of the four-momenta p and Q 
and polarizations and v a have to be made. The functions fj,j±i(q 2 ) and fj,j±i{q 2 ) for 
the P**-states are as follows: 

W* 2 ) = 2 2 2 ff* (19) 

m%-p 2 - q 2 + 2ajj ±1 ^p 2 q 2 

/W? 2 ) = —~ 2 ^ t=- (20) 

q 2 + mr D * -m 2 + 2a J)J ± 1 m D , yf q 1 

5 Results 

Using the technique described in the previous sections, we present the final expressions 
for matrix elements with different intermediate resonances. The total matrix element 
squared is as follows: 

\M\ 2 = \M 6 + M po + M p(1450) + M 6l(1235) + M b2 + M p2 + M P3 + M Dl + M D{ + M D * + 
+ M- +M D2 +M D ,+M 3 - | 2 . (21) 

x 3/2 2 °5/2 

Here, M$ presents the non- resonant contributions to the matrix element. The amplitudes 
Md 1 and Mr>> are as follows: 

M Dl = sin^ \ 3u = 1/2 > + cos^x \j u = 3/2 > , 

M D[ = cos^i \j u = 1/2 > - sintfi e i,?2 |j u = 3/2 > , (22) 

where $i and $2 are mixing angles and similar expressions can be used for 2~-states. 
The resonant matrix element can be presented as follows: 

M Rj = ^V^ ^'^y^ £ C Ll C L2 F L M 2 )hM 2 )V LlL2 A LlL2 . (23) 

L1L2 

Here, Li(L 2 ) is the angular orbital momentum in the B°(R) rest frame; Cx 13 Cx 2 are 
relative amplitudes defined above, Vl 1 l 2 is the expression for the momentum dependence; 
At, x L-2. is the expression for the angular dependence. The expressions Vl 1 l 2 and Al x l 2 are 
combined in Table [1] for different intermediate states. The notations c a = cos a and 
s a = sin a are used. The functions fj,j±i(q 2 ) and /j,j±i(<? 2 ) used in Tabled] are defined 
by (lT3"j) and ffTSj) for the co>7r-resonances and by <HM and (1201) for the P**-resonances. 
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-l^SgS^SpS^ - SgCtpSpC^C 2 - - 1/2) 


T~ P 

5/2 


F 


— 7r m B/2,3(g 2 )Qp 


-3/2c e c^(c? - 1/3) + - I/Isqc^c^- 
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SftsoSfCps^ + sgs^spc^c^ - sgc^s^s^ 


F 


F 


4 ? mB/2, 3 (9 2 )/ 2 , 3 (9 2 )Qp 


9/Ac0C^(c^ - 1/3) + 3/4c0S0C^s 2 £ + S/Asgc^s^- 
-l/ZseSfSps^ct - sec^spc^c 2 - - 1/2) 


3 5/2 D 


F 


-j=m B Q 3 p 2 

v<? 2 


seCsSas^ctsl - 4/15s0Cri,s / 3S V ,C£+ 

+ l/3s 6 >S4,S /9 C, / ,(c| - 1/5) - CgSpS^(cj - 1/5) 


F 


F 




(s^s^s^ + s e c^spc^){(^ - 1/5) - 2s0S ( £S /9 s v ,C£s| 


G 


t 


i ™2 x /'„2\r~»3„2 


sec^s^s^sj - 4/15s0C0S / 3S^,C£+ 

+ l/3s e S0S /3 C^(c| - 1/5) + 4/3c S/3S^ s ?(cf - 1/5) 



Table 1. Summary of momentum and angular distributions 
for different intermediate states, which are described in this 
paper. 



6 Decay chain simulation 



To demonstrate the angular distributions for each intermediate resonance in the D*ujix 
final state, we generate 2 x 10 6 B° — > D* + uon~ events according to the phase space 
distribution using the qq98 program package p2] . For a further study we fill profile angular 
spectra with the appropriate weight density functions for each resonant hypothesis, which 
have been obtained above. 

A description of each vertex includes transition form factors. Since it is not yet 
possible to obtain these form factors from rigorous theoretical calculations, we rely on 
the simple phenomenological Blatt-Weisskopf model [TSJ [TH]. For L > this simple 
form factor suppresses growth of the matrix element with final particle momentum. The 
Blatt-Weisskopf functions Bl(x) are chosen as follows: 

Bl{x) ~ x^\h L {x)\ ' (24) 

where 

h L (x) = Z! e *(«-¥) £(-i r ^-±^(2u;r (25) 
x nUL — n)\ 

n=0 v ' 

is a spherical Hankel function, x = kr, x = kor, k, k are the magnitudes of the daughter 
particle three-momentum in the mother particle rest frame for the case when the resonance 
four-momentum squared is equal to q 2 and m 2 R , respectively, and r = 1.6 GeV -1 is a 
hadron scale. According to our normalization, these functions are equal to one, when 
y^q 2 = tur. Another common normalization gives Bl(x) = 1 for x — 1. The Blatt- 
Weisskopf functions corresponding to L discussed here are given below for convenience: 



B Q (x) = 1, 
B x {x) = 



1 + x 2 



1 + x 



_0 

.2 : 



_ (xl 3) 2 + 9xq 
B2[X) ~ 3)2 + 9*2' 



B 3 (x) 



jx 2 {x 2 - 15) 2 + 9(2xg - 5) 2 
x 2 (x 2 -15) 2 + 9(2x 2 -5) 2 ' 



- lo.r 2 • l()o) 2 • 2o,- 2 i2.r 2 -21) 2 
" W fx 4 - 45x 2 + 105) 2 + 25x 2 (2x 2 - 21) 2 ' 1 J 
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Figure 5. Simulated angular distributions for the w7r-resonances. The figures al), bl), cl), dl), 
el) correspond to the J p = 0~ {Pq) intermediate state; a2), b2), c2), d2), e2) — J p = 1~ 
(p(1450) )-state; a3), b3), c3), d3), e3) — J p = 1+ (6i(1235) _ )-state. 
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Figure 6. Simulated angular distributions for the 0J7r-resonances. The figures al), bl), cl), dl), 
el) correspond to the J p = 2 + (b^) intermediate state; a2), b2), c2), d2), e2) — 2~ (p^)-state; 
a3), b3), c3), d3), e3) — 3" (p 3 (1690)-)-state. 
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V V V 

el) e2) e3) 

Figure 7. Simulated angular distributions for the P-wave Z)**-resonances. The figures al), bl), 
cl), dl), el) correspond to the Jj^ = 1^ 2 narrow state; a2), b2), c2), d2), e2) — = ly 2 

broad state; a3), b3), c3), d3), e3) — = 2+, 2 narrow state. 
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Figure 8. Simulated angular distributions for the D-wave Z)**-resonances. The figures al), bl), 



cl), dl), el) correspond to the Jj u 

state; a3), b3), c3), d3), e3) 
narrow state. 



l 3 y 2 broad state; a2), b2), c2), d2), e2) 



J u 



2 3 y 2 broad state; a4), b4), c4), d4), e4) 
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If we consider one angular variable only, the distributions can be the same for different 
resonant hypotheses. Efficient separation between resonances is possible, when all angu- 
lar variables are taken into account. This statement is demonstrated in Figs. [5] and |6] for 
the uii- states and in Figs. [7] and M for the _D*7r-states. As mentioned above, J p = 1 + 
P-wave and J p = 2~ D-wave states are a mixture of pure states. However, for demon- 
stration purposes we consider and show angular distributions for pure states. Moreover, 
we use a simple relativistic quark model of mesons to estimate constant ratios Cj^i/Cj 
and Cj + i/Cj in pip , which are responsible for relative contributions of amplitudes with 
different orbital momenta in the total matrix element [20]. The constant ratios for all 
discussed states are chosen roughly as follows: Cj^i/Cj = 3/2, Cj + i/Cj = 2. 

For Dalitz plot analysis, interference between resonances should be taken into account. 
For a one-dimensional distribution, an interference term for resonances, which decay to 
the same final state, can cancel out after integration over other variables. However, 
in a real experiment such cancellation can disappear due to the nonuniform detection 
efficiency, so that a finite interference term can be observed. For resonances, which decay 
to the different final states uir and D*tt, the interference term cannot be neglected. For 
demonstration purpose we show distributions between 6 1 (1235) _ and pure as well as 
p(1450)~ and pure D®. However, there is possible interference between the resonant and 
non-resonant structures. 

For simulation we use BW functions for 6i(1235)", p(1450)~, L>° and D?. Thus, the 
g 2 -dependent widths have to be obtained. For a g 2 -dependent width of the b\ we consider 



19 



the dominant decay to the ujtt [2]: 



i 6iw J - -7=^ 2 — ; — ?9 / 9 \ ^'J 

vV m^ + bofl^mij Po 



where po is the magnitude of the u momentum in the resonance rest frame, when q 2 = m 2 r 
Here, we use the fact that the experimental ratio of the amplitudes with L = 2 and L = 
is about 0.3 [2] and thus the constant bn ~ 53. For a g 2 -dependent width of the pure Di 
and pure D[ we consider the decay to the D*tt: 

TzW^ 2 ) = ^^ 2 )£r Dl , D ,, (28) 
y q 

where Qo is the magnitude of the D* momentum in the resonance rest frame, when 
q 2 = m 2 D D , . For a g 2 -dependent width of the p(1450) we consider its decays into the 
a!(1260)7r and utt modes: 

n , m p(1450) m* {im F 2 (q 2 ) + b D f 2 2 (g 2 )F 2 (q 2 ) k(ai) 

i (9(1450) - I-L-OJ -j= " — ~ " lp(1450) + 



m J(i450) +bDfl 2 (m 2 piU50) ) k 0(ai) 
+ a ^ F 2 (g 2 )H!r p(1450) , (29) 

"V(1450) Po 

where the parameter a = 2/5, when ^J'q 2 > m ai + m,r- [21J, and a — 1, when y^g 2 < 
m ai + m^-, k( ai ) is the momentum of the a x in the p(1450) rest frame, k ( ai ) is the same 
momentum, when i/g 2 = ^(1450) • Here, we use the fact that the experimental ratio of 
the amplitudes with L = 2 and L = in the ai(1260) — > pn decay is about —0.06 [2] and 
accept the same value for the p(1450) — > a!(1260)7r decay. Thus, we can estimate the 
constant be ~ 182. 

Obviously, it is impossible to analyse spectra without knowledge of the relative phases 
in the amplitude. Thus, in Fig. [9] we show some typical distributions for different relative 
phases Aip, such as 0, tt/2, tt, 37r/2 and the distribution without interference. The relative 
constant amplitudes between resonant matrix elements squared are chosen of one order 
of magnitude for the 61 (1235) ~ and D\ for simplicity and one order of magnitude smaller 
for the Di than for the p(1450)~ according to experiment [12]. Although small, the 
interference effects are not negligible. 
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bl)p'D' 1 b2) p'D[ 



Figure 9. Demonstrative interference distributions. The figures al), a2) correspond to the 
distributions over angles cos£i and cos^2 for interference between the 6i(1235) _ and pure D®; 
bl), b2) correspond to the distributions over angles (pi and <p2 for interference between the 
p(1450)~ and pure D®. The subscripts 1 and 2 correspond to the ujtt- and Z)**-resonances, 
respectively. 



7 Conclusion 

We have described a model of the B° — > D* + um~ decay, in which a total amplitude is a 
sum of contributions of different intermediate states. In our study we consider different 
resonant contributions to the matrix element, such as light 0J7r-hadrons with the spin- 
parities of J p = 0~, 1*, 2 1 * 1 , 3~, and heavy-light hadrons, which are excitations of the cu- 
states in P- and D-waves. All resonances are described by the relativistic Breit-Wigner 
factors. The resonant matrix elements are parameterized in the angular basis, which is 
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convenient for the experimental Dalitz plot analysis and is natural from the physical point 
of view. 

Monte-Carlo simulation based on the obtained expressions has been performed. The 
angular distributions obtained for the listed above intermediate states and their interfer- 
ence effects are demonstrated. 
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A Appendix 

In this section we present the phase integral W{p 2 ) at the decay rate defined by (ITDl . 
The integral 



W{p 2 



7T 



(2m+) 2 



"^4- max 

dM 2 I dM 2 



a 3 n + ^2 



QpTTTT 



(Al) 



is a standard phase space factor for w-decay [22J. Here, the Kibble determinant A{p, P + , P ), 
which zeros determine the phase-space boundary, is presented as follows: 



A(p,P + ,P ) 



p 2 pP + pP 
pP + m 2 P + P 
pP P+P m 2 



where m and mo are the charged and neutral pions masses, respectively; the range limits 
of Mi are 



Ml 



M, 



(E- 



(+-) , E (+-)\2 



(p| + - } + E t 



(+-)\2 
) 



where 



E. 



(+-) 



E, 



(+-) 



{\lE^~ )2 -m 2 



Pi + ~ )2 -m 2 



M 2 - 2m 2 

2~Mo ' 
p 2 - M 2 - m 2 

2Mn 



A) - m 0) > 



E, 



(+-) 2 ^2\2 







m 2 r, (A2) 



(A3) 



are the energies of ir + and 7r° from u decay in the 7r + 7r rest frame; 

. 9x 1 + (m^r) 2 



9u/m(p 



Quipir \ m u 



1 + (\/p 2 r) 2 



(A4) 
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is the form factor which restricts too fast growth of the width T^p 2 ) with p 2 , so that 
To;(p 2 ) — > const as p 2 — > oo (here r is a hadron scale) [23j; the quantity <7 W/OT (?™ 2 ) — 
16 GeV -1 [2U [251 [26]; the and g pnn amplitudes are assumed to be real constants and, 
thus, a 37r = (0.01 ± 0.23 ± 0.25) x 10" 5 MeV" 2 [22] and g p7Tn ~ 6 



D p ±,o (Mj_ )0 ) = M 2 >0 - m 2 ± , + m, ± , r p± ,o (M 2 10 ) , ( A5) 

where 

M 2 = p 2 -Ml- M 2 + m 2 + 2m 2 , 



3 



r ^' o(M± ' o) = (^K^y 1 r ^ ' (A6) 

k{ is an absolute value of pion momentum in the 7r + 7r~ rest frame for i = 0, in the 7r + 7r° 
rest frame for i = + and in the 7t~tt° rest frame for z = — ; 

Z(M ±:0 ) = 1 - zsi$(M±, , V?) (A7) 

is the factor taking into account the interaction of the p and 7r mesons in the final u decay 
state, where the parameter si = 1 ± 0.2 corresponds to the prediction of [27], where the 
specific form of the $(M ±)0 , a/p 2 ") function can be found. The couplings g upw and g p7T7T are 
the same for p ±,Q because of isotopic invariance. As emphasized in the text, the angle £ 
is related to the intermediate resonance mass by a simple expression. 

Finally, let us give here the p 2 -dependent width of the w-meson taking into account 
the 37r and 7r°7 modes [2]: 

r (n 2 \ - w{p2) b r i < (p 2 -<) 3 9L,(p 2 ) ( 

^ [P) - W(mlf^ L " + p 2 (ml - mlfgl^mlf^"- { ^ ] 
A form factor 5WyGo 2 ) has a similar form ( 1A4I) [22] and gunyiWu) — 0.7 GeV -1 [26J. 
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